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Abstract. In materials that do not allow birefringcnt phase- matching or 
periodic poling we propose to use waveguides to exploit the tensor structure of the 
second order nonlinearity for quasi-phase matching of nonlinear interactions. In 
particular, we concentrate on curved waveguides in which the interplay between 
the propagation direction, electric field polarizations and the nonlinearity can 
change the strength and sign of the nonlinear interaction periodically to achieve 
quasi-phase matching. 

Nonlinear optics enables many technologies in classical and quantum communi- 
cation. The best known application is second-harmonic generation (SHG) pQ, the 
frequency doubling of an electromagnetic wave. Quantum mechanically the reverse 
process can happen as well and a single photon can split into two photons under 
energy conservation. This process is called spontaneous parametric down-conversion 
(SPDC) and has been successfully used for generating entangled photon pair states 0, 
or heralded single photons 3,4. SPDC is the workhorse for quantum communication 
[3 [6] and optical quantum information processing [7] . 

The second order term in the perturbative expansion of the electric polarization 
caused by a driving field is responsible for SHG, sum- and difference frequency 
generation (SFG and DFG), and SPDC. In all of these processes waves with vastly 
different frequencies interact. Dispersion causes the generated waves to de-synchronizc 
from the driving field. This is detrimental to the growth of these new frequencies, 
and one has to take special measures to avoid this and achieve phase-matching. In 
birefringent materials one can phase-match by choosing the polarizations, crystal cut, 
and crystal temperature so that the birefringence effectively cancels the dispersion. 
Furthermore, by structuring the material into waveguides with suitable modes [5], by 
introducing form birefringence [9j 110) one can achieve matching effective refractive 
indices. 

All other techniques employ some form of quasi-phase matching (QPM) by 
judiciously introducing periodicity to rephase the waves |11) . Photonic crystals offer 
the highest flexibility by introducing periodicity both in the linear and nonlinear 
properties of the material [12] . Traditional QPM methods only periodically reverse the 
sign of the nonlinearity. The most common implementation is by periodic poling of a 
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ferroelectric. The period is chosen so as to reverse the sign just when the accumulated 
phase difference begins to drive energy back into the exciting field. The sign reversal 
effectively resets the phase and allows the conversion to continue. Waveguides have 
the added benefit that the light stays contained and other optical elements can be 
integrated into the structure. In this article, we introduce another kind of quasi- 
phase-matching that employs waveguides. We call this method directional quasi- 
phase-matching (DQPM). 

DQPM makes use of the fact that the strength of the nonlinear interaction is 
governed by a tensor x\jk ~ or dij m ^ ne usua -l contracted notation |13) . Therefore, the 
coupling strength between three waves depends on the orientation of the three fields 
with respect to the crystal and each other. If we assume that the field propagation 
direction in a waveguide adiabatically follows the direction of that waveguide and the 
field polarizations stay approximately transverse to the propagation direction, then 
the coupling can be modulated by changing the waveguide direction. This modulation 
can affect both the sign and magnitude of the coupling. To our knowledge, apart from 
some earlier work on bulk materials |14) this effect was only noted very recently in 
connection with nonlinear optics in disk and ring resonators |15[ 1161 117] . 

In what follows, we demonstrate the conditions for DQPM as they apply to SHG 
and the related phenomena in curved waveguides. We discuss various curve shapes in 
relation to the expected amplification of SHG, their bend loss and modal behavior. 
Our analysis also applies to other material systems - in particular other compounds 
that crystallize in the zincblende structure - but presently we focus on the Al x Gai- x As 
material system. Al x Gai- x As is a common optoelectronic material that is transparent 
for telecommunications frequencies and has a high optical nonlinearity (100 pm/V). 
Additionally, semiconductor lasers and other active components are manufactured in 
this system, making it an ideal choice for further integration. 

In our analysis, we begin with the standard treatment of SHG in straight 
waveguides and modify the equations to account for directional changes. We neglect 
losses and use the slowly varying amplitude approximation. We further assume that 
the fields are mainly transverse to the direction of propagation. For bulk Al x Gai_ x As 
which belongs to the 43m point group, the only non-zero elements of the second order 
nonlinear tensor are ^14=^25=^36, where, as usual the labels are 1,2,3 = x,y,z and 
1,2,3,4,5,6 = xx,yy, zz,yz,xz,xy in the first and second indices, respectively |13j . 
This means that the two fundamental field modes and the second harmonic are only 
coupled if there are nonzero field amplitudes along every one of the crystallographic 
axes x, y, and zX 

To develop a visual sense of DQPM, consider two co-ordinate systems. One is 
static and oriented with respect to the cubic symmetry of the AlGaAs lattice, the other 
co-moves with the direction of propagation of light in the waveguide. For Al x Gai_ x As, 
epitaxial growth is typically done on (100) surfaces (y — z plane), therefore a waveguide 
is confined to lie in this plane. We denote the waveguide propagation direction tangent 
vector as ds as shown in figure Q] If ds is parallel to z, then all guided fields will 
normally be polarized along y or x. In this case the effective nonlinearity will be 
zero, and there won't be any nonlinear conversion, because none of the fields have a 
z component. On the other hand, if the waveguide is directed at ±45° to the y-axis, 

I In tightly confining Bragg reflection waveguide 10 samples provided by A. Helmy we have recently 
observed experimentally that other interactions (e.g. all polarizations equal) are possible that 
wouldn't be allowed in the bulk material. Such interactions may occur because the symmetry is 
broken by the layer structure of the respective waveguides. 
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Figure 1. Crystal coordinates, waveguide tangent vector ds and its angle a with 
the y axis. 



say ds oc y + z, then its two transverse polarizations are parallel to x and parallel to 
y + z. In this situation the nonlinear process can, at least in principle, occur. 

We will now use the dynamical equations for the amplitudes of plane waves in 
second harmonic generation, assuming that these amplitudes vary slowly compared to 
the wavelength. The plane wave approach is justified as long as the waveguide modes 
that are involved are mostly transverse, however, the coupling strengths will be off 
because the transverse mode structure is being neglected here. If s is the arc-length 
parameter of the waveguide curve, the change in the amplitude of the second harmonic 
field E2(s) is given by 

±E 2 (s) = -iO* (1) 
as n 2 c 1 

The subscript indices 1 and 2 represent the fundamental and second harmonic fields 
respectively. These fields are defined by their respective amplitudes (E), frequencies 
(<jj), and wavevectors (k), and the phase mismatch Ak = k 2 — 2ki > 0. n 2 is the 
effective refractive index for the second harmonic field. 

We consider the generation of vertically polarized (x) second harmonic light (E 2 ) 
by injecting horizontally (in the y — z plane) polarized fundamental light (Ex) into the 
guide. Therefore we evaluate the effective nonlinearity d c s as 

defi = d>ijka 2 iaijaik (2) 

ijk 

= di4,a 2x ax y a\z = d\4,ax y a\ z , (3) 

where a\ and a 2 are unit vectors in the polarization directions of the fields E± and E 2 
respectively. Further a 2x = 1 because the polarization of the second harmonic field is 
always along x. 

We now consider arbitrarily directed waveguide structures and reinterpret s as 
the effective path length along the waveguide. In the simplest approximation s will be 
measured along the midline of the guide. For a more realistic treatment one would use 
the corrected propagation constants for the modes of a bent waveguide, or short of a 
full two- or three-dimensional simulation, use the path taken by the field maximum. 
We retain the adiabatic slowly varying amplitude approximation for waves that could 
either be plane waves or confined modes with quasi-transverse polarizations. Also 
we assume that the magnitude of the phase mismatch between the fundamental and 
second harmonic waves remains constant (Afc ^ Afc(s)). 

We define the local direction of the waveguide by its angle a(s) with the y-axis, 
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Figure 2. First order directional quasi-phase matching. In this highly idealized 
structure we assume no phase shifts, loss, or out of plane field rotations occurring 
at the sharp 90° bends. The change in second harmonic amplitude (right side of 
(0) switches from positive to negative because of the sin(2a(s)) term, depending 
on the quadrant of a. For first order QPM each segment have to be exactly L c 
long, so that the sign change occurs precisely when, due to the dispersion, the 
phase fronts between the pump and SHG fields have acquired a it phase between 
each other. 



resulting in an expression for the in-plane polarization vector ai that is given by 

I ° \ 

a i( s ) = — sina(s) . (4) 
\ cosa(s) J 

Making the substitutions for d we get 

^ =i ^ sin[ 2a( s )]^ e< (A fe ), (5) 
as n^c 

The growth of E 2 is now clearly dependent on both the distance traveled and 
the orientation of the waveguide, a(s). It is zero for a waveguide whose propagation 
direction is along any of the crystal axes and can have positive or negative sign. 

We will assume that the exciting field E\ has constant amplitude, i.e. is not 
depleted, and that ^(s = 0) = 0. For straight guides with sin2a(s) = const. > and 
Afc(s) ^ 0, (0 describes the familiar behavior for non-phase-matched SHG. That is, 
E2 oscillates with a period of 2L C where L c = Tt/Ak, between zero and it's maximum 
value along the guide. 

DQPM can thus mimic the behavior of the usual periodic poling of a ferroelectric 
nonlinear crystal. In a "zig-zag" shaped guide with segments of length L c that 
alternate between the a = tt/4 and a = — 7r/4 directions the nonlinear polarization 
switches its sign as in periodic poling. Figure [2] shows the layout of the structure on 
the wafer. 

Clearly such sharp bends would be difficult, if not impossible to make, given that 
L c for the doubling of light at 1600 nm wavelength to 800 nm in bulk Alo.4Gao.6As 
is just under two microns. Instead of bends it would be better to use total internal 
reflection mirrors, etched into the sample |18j . with an appropriate treatment of the 
phase change upon reflection. In any case, the losses would likely be very high, 
therefore we will now consider perturbations of first order DQPM as well as higher 
order phase matching to allow for milder bending. The fact that one has more control 
over the effective nonlinearity in DQPM than in poling reflects in the fact that even 
order QPM is possible here. As an example, we will look at a second order QPM 
scheme, with the main advantage that it only requires 45° bends. 

Figures [3] and H] depict second order and perturbative DQPM respectively. These 
designs are intuitive because their subsections are in units or half units of L c . 
This makes it relatively easy to keep track of the phase accrual between the pump 
and second harmonic fields. The second order structure turns off the non-linearity 
whenever it's propagation axis lies along one of the crystallographic axis. Compared 
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Figure 3. This structure is an example of an even order quasi phase matching 
scheme - something that is unavailable with periodic poling methods. This is 
because poling cannot turn off the non-linearity, only reverse it. It also allows 
somewhat milder bending. However, similar growth of SHG to the first order 
scheme will require twice the distance. 




Figure 4. These two DQPM structures have no kinks and allow studying a 
pcrturbative regime with only small deviations from the straight line. Therefore 
they will in practice allow balancing bend losses and nonlinear gain. (Top:) One 
period of this structure is composed of four circular arc segments that are each 
L c /2 long. For large radii they will only attain small angles off the horizontal. 
(Bottom:) This structure is a cosine-shaped curve with a path length of 2L C . 
Because it has a continuously changing curvature it avoids the mode overlap losses 
that can happen in waveguides with a circular structure. However, at radii that 
are large with respect to the waveguide width the two structures become almost 
indistinguishable. 

to first order DQPM, second order DQPM requires twice the distance to generate a 
similar intensity in the second harmonic. 

The perturbative structure consists of two oppositely directed S-bends, each of 
length L c , placed in series with each other. Each S-bend consists of two circular 
segments of length L c /2. This produces a structure with smoothly varying angle 
a(s) and a curvature of constant magnitude but alternating sign.^ Over the first S- 
bend, the transverse electric field direction stays in the quadrant where sin(2a(s)) is 
positive. Over the second S-bend, the electric field direction stays in the quadrant 
where sin(2a(s)) is negative. As in first order DQPM, this coincides with respective 
positive and negative contributions from the accumulated phase between the pump 
and second harmonic fields. Thus, the net result is always a positive contribution to 
the overall field amplitude of the second harmonic. The mild bending doesn't produce 
nearly the amount of polarization response in the medium and thus the conversion 
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Figure 5. Growth of the second harmonic field amplitude vs. the path length 
along the guide over two periods for first and second order DQPM, for the circular 
perturbative DQPM with R = 5L C , and for the cosine-shaped perturbative DQPM 
with R > 5L C . 



efficiency to the second harmonic is much less. 

Depending on the precise modal structure of the waveguide a sudden change in 
curvature could be inducing more loss than an adiabatic one. Therefore we investigated 
cosine-shaped curves as alternatives to the segmented circular S-bend (see figure [4]). 
In a cosine curve the curvature changes smoothly along the curve and thus mode 
overlap losses can be avoided. For large curvature radii and small lateral displacement 
the circular and cosine shaped S-bends become very similar both in shape and in 
performance. 

Each design is scalable by virtue of the 2n relative phase between fundamental 
and second harmonic introduced by one unit cell and every period yields the same 
incremental amplitude in the undepleted pump approximation. This will lead to an 
overall quadratic growth of the second harmonic intensity. Figure [5] plots the value 
of the integral (ignoring the constant prefactors) in ([5]) for two periods of each of 
the four structures. While (O can be integrated analytically for all four structures 
with elliptic integrals for the cosine shaped curves, we used numeric integration to 
obtain these results, because of the greater flexibility. It is obvious that the more 
realistic, perturbative designs fall quite short of the maximum possible conversion 
efficiency for reasonable curvatures. With strong index contrast as, for example, in 

§ A special case is a full circle, as it would occur in a ring resonator. In this case we would first have 
to solve for the ring resonance condition at the two frequencies and then check the coupling between 
the two modes, which is modulated simisoidally around the ring with a 90° period. Ignoring the ring 
resonator modes for a moment, we can integrate {5} around a ring with a circumference 27rr = NL C 
that is an integer multiple of the nonlinear coherence length. We find that for N ^ 4 the field at 
2io will always be identically zero after two loops. When N = 4 a quasi-phase matching condition 
is achieved and the field will continue to grow. It may appear surprising at first that no higher- 
order QPM is possible, but if we realize that the ring yields a sinusoidal modulation of the effective 
nonlinearity, it is clear that there are no higher Fourier coefficients available for QPM. Therefore any 
ring-resonator QPM design seeks to use mode pairs with large L c , in order to keep the ring radius 
reasonably big |16l 117] 
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Figure 6. SHG field amplitude after one period of the circular perturbative 
structure as a function of the curvature radius. 

deeply etched ridge guides, these curvatures still seem feasible and it will remain to 
be shown by experiment how the nonlinear gain compares to the losses. Figure [5] 
shows the amplitude gain of one period of the circular structure as a function of the 
curvature radius R. 

The ideal zig-zag of first order DQPM involves sudden changes in waveguide 
direction, and possibly violates some of the assumptions that underpin the basic 
analysis. Whether mediated by sharp bends or reflections, the design would almost 
certainly attenuate the field amplitudes after each successive direction change and 
possibly add a certain phase shift in the reflections.-^ 

One of our simplifying assumptions is that the local TE polarization follows the 
normal to the waveguide centerline. It is clear that this cannot be true in the real 
situation. In order to get a feeling for the actual field path, we have performed 
two-dimensional Finitc-Differcncc-Timc-Domain (FDTD) simulations of the circular 
waveguide structure shown in the top part of figure HI using a commercial solver 
(Lumerical, FDTD Solutions 6.0). Figure [7] shows the power distribution, local 
polarization and the polarization angle along the field maximum. While we find 
that the field maximum takes a path that is quite different from the centerline, the 
polarization stays almost perfectly transverse, except near the edges of the guide, 
where it is comparably weak. This means that by integrating ([5]) along the actual 
field path one could design the ideal curvature and length to optimize the conversion 
efficiency. In practice, it would certainly be easier to use FDTD solvers that can 
treat the second order nonlinearity perturbatively while propagating the field and 
thus provide a complete simulation of the device. It is remarkable that the maximum 

|| For total internal reflection, this phase is precisely it, and for three waves, each wave would phase 
shift for a total change in the argument of the exponential of AfcL c + 3n = 4-7T. This would mean 
that the phase has essentially already been reset by total internal reflection. An additional ir phase 
from the ensuing direction change in the guide would actually completely negate the SHG process 
This effect has been suggested as a possible solution for QPM but the idea of potentially different 
coupling on the direction change was not mentioned. For typical single mode waveguides in AlGaAs, 
L c is quite small, and to do effective phase matching in this manner would be quite difficult. One 
would need to either effect nearly 180° reflections or very shallow reflections in order to ignore the 
effects of DQPM. Further it is unclear that the physics of total internal reflection would apply to 
standard waveguide modes. 

In connection with strongly sloping sidewalls sharp bends can also induce facilitate the polarization 
rotations | 19l 1201 121| . With this one could conceivable also achieve some kind of QPM if the pump 
field could be made to rotate from in plane to out of plane with the correct periodicity. 
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Figure 7. Linear FDTD simulation of the circular perturbative DQPM structure. 
Top: False color plot of the intensity of TE (in-plane) polarized electric field at a 
vacuum wavelength of 1550 nm propagating through one period of the structure. 
The waveguide has refractive index 3.0, a width of 0.5 (im that and is surrounded 
by air. Here, the field is being injected in the fundamental waveguide mode 
from the left. It shows negligible loss but considerable field intensity outside the 
waveguide edges (thick white lines). Bottom left: Quiver plot of the local field 
polarization. The fact that this is always orthogonal to the local propagation 
direction can be visually verified along the path of the field maximum (thick 
green line). Bottom right: Local propagation direction (solid blue line) and the 
respective local polarization direction (—90°, dashed magenta line) at the field 
maximum as a function of the horizontal coordinate along the waveguide. For 
comparison, the angle of the tangent to the waveguide ccntcrlinc is shown with 
the dash-dotted green line. 



curvature of the field path is quite a bit stronger than the waveguide's curvature. This 
means that with an adapted design we could expect a higher conversion efficiency than 
in our simple case. 

The tight bends required by our proposed structures may seem very challenging, 
but a lot of theoretical and experimental research has been performed on this topic. 
Very tight bends [H[23], including 90° bends [24] and more adiabatic bends [251126] . 
have been investigated. Important observations are: 1) the effective index decreases 
with decreasing bend radius; 2) wider guides exhibit more drastic modal changes in 
bends than do narrower guides; 3) Strong radial confinement is necessary to inhibit 
radiative loss in a bend; 4) high radial symmetry is necessary to reduce intermodal 
mixing. In essence coupling to radiative (non-guided) modes as well as intermodal 
coupling is inevitable in any waveguide structure that incorporates bending. 

We have assumed that the mode profiles of the three interacting fields stay the 
same and that the ideal structure would enforce single mode behavior everywhere 
along it's length. In order to avoid intermodal mixing, excellent control over the 
spatial profile of the waveguide is desirable, which will typically result in very high 
aspect ratios. While this kind of structure can be fabricated, it requires state of the 
art manufacturing techniques that are extremely dependent on the proper chemistry 
during the fabrication process. Unfortunately, these processes also often have an 
unwanted side effect - sidewall roughness. This will be the main source of loss due 
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Figure 8. A DQPM design featuring the perturbative scheme. Here there are 
13 unit cells repeated, and the perturbation is quite large so that bend losses 
are expected to be significant. We find that light is being guided through the 
structure but we only have very preliminary loss estimates. 



to the high index contrast and tight mode confinement for waveguides that have been 
etched through the core. 

Success of these designs will thus depend on how the conversion efficiency- 
compares to the various loss mechanisms. As there are in principle an infinite number 
of variations that might work, the actual design must incorporate fabrication methods 
as an input variable. A first fabrication try for a DQPM structure is shown in figure[5] 

Modeling must now be performed to best predict design behavior based on the 
available fabrication capabilities. In the future, we hope to make considerations for 
such related issues like etching and regrowth processes and their effects on insertion 
and bend loss. Numerical simulations are underway to predict optimal mode matching 
between waveguides of different curvature and optimal aspect ratios. Lastly, the basic 
theory may need some modification to allow for modal conversion along the waveguide 
structure. We believe that directional quasi-phase matching is an interesting and 
very flexible route towards exploiting second-order nonlinearities in many important 
semiconductor materials. These materials are notoriously hard to phase-match. Our 
concepts can be generalized to non-circular resonators, such as stadium resonators, in 
order to exploit higher-order phase-matching together with the increased field strength 
afforded by a cavity. 
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